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Abstract 

We extend the hfecycle model (LCM) of consumption over a random horizon (a.k.a. the Yaari model) to a world 
in which (i.) the force of mortality obeys a diffusion process as opposed to being deterministic, and (ii.) a consumer 
can adapt their consumption strategy to new information about their mortality rate (a.k.a. health status) as it becomes 
available. In particular, we derive the optimal consumption rate and focus on the impact of mortality rate uncertainty 
vs. simple lifetime uncertainty - assuming the actuarial survival curves are initially identical - in the retirement phase 
where this risk plays a greater role. 

In addition to deriving and numerically solving the PDE for the optimal consumption rate, our main general result 
is that when utility preferences are logarithmic the initial consumption rates are identical. But, in a CRRA framework 
in which the coefficient of relative risk aversion is greater (smaller) than one, the consumption rate is higher (lower) 
and a stochastic force of mortality does make a difference. 

That said, numerical experiments indicate that even for non-logarithmic preferences, the stochastic mortality ef- 
fect is relatively minor from the individual's perspective. Our results should be relevant to researchers interested in 
calibrating the lifecycle model as well as those who provide normative guidance (a.k.a. financial advice) to retirees. 

Keywords: lifecycle consumption, stochastic mortality, survival curve matching, JEL codes: E21/G22, Subject 
Category: ??, Insurance Branch Category: ?? 



1. Introduction and Motivation 

The lifecycle model (LCM) o f saving s and c onsump- 
tion - originally postulated by 'pishe 7 ('1930') and re- 
fined by Modigli ani and Brumb erg (1954), Modigliani 
(fl986l) - is at the core of most multi period asset pricing 
and allocation models, as well as the foundation of mi- 
croeconomic consu mer behavior. T he orig inal formu- 
lation - for example iRamsev ( 1928 ) and lPhelpsI ( 1962 ) 
- assumed a deterministic horizon. But, in a seminal 
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cont ribution, the LCM was extended by lYaari (Il96i 
Il965i) to a st ochastic hfetim e , which eventua lly led to 
the models of lMertonI (Il97lh . iRichardI d 1975b and hun- 
dreds of subsequent papers on asset allocation over the 
human lifecycle. 

The conceptual under-pinning of the LCM is the in- 
tuitive notion of consumption smoothing whereby (ra- 
tional) individuals seek to minimize disruptions to their 
standard of living over their entire life. They plan a con- 
sumption profile that is continuous, equating marginal 
utility at all points, based on the assumption of a con- 
cave utility fu nction. See the r ecent (and very accessi- 
ble) article bv iKotlikoffi (120081) in which this concept is 
explained in a non-technical way. 



Once agai n, until the seminal contribution by lYaari 
(Il964lll965h . the LCM was employed by economists in 
an idealized world in which death occurred with prob- 
ability one at some terminal horizon. Menahem Yaari 
introduced lifetime uncertainty into the lifecycle model, 
in addition to - his more widely known contribution of 
- introducing actuarial notes and annuities into optimal 
consumption theory. 
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In the expressions (and theo rems) he der ived for the 



optimal consumption function. I Yaaril (Il965h assumed a 
very general force of mortality for the remaining hfe- 
time random variable, without specifying a particular 
law. His results would obviously include a constant 
force of mortality (i.e. exponential remaining lifetime) 
as well as Gompertz-Makeham (GM) mortality, and 
other commonly formulated approximations. Yaari pro- 
vided a rigorous foundation for Irving Fisher's claim 
that lifetime uncertainty effectively increases consump- 
tion impatience and is akin to behavior under higher 
subjective discount rates. Mathematically, the mortal- 
ity rate was added to the subjective discount rate. 

That said, most of the empirical or prescriptive pa- 
pers in the LCM literature have not gone beyond assum- 
ing the GM law - or some related deterministic func- 
tion - for calibration purposes. In other words, mortal- 
ity is just a substitute for subjective discount rates. In 
fact, one is hard-pressed to differentiate high levels of 
longevity and mortality risk aversion from weak pref- 
erences for consumption today vs. the future, i.e. pa- 
tience. Some have labeled this risk neutrality with re- 
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Feieenbaumj (2008|). or the recent work bv 'Lachance 
( 2012) - all assume a deterministic force of mortality. 



Indeed, some economists continue (surprisingly) to 
ignore mortahty a lltogether, for example the recent re- 

Q). Perhaps this is 



view by lAttanasio and Weber 



because when the force of mortality is deterministic, it 
can be added to the subjective discount rate without any 
impact on the mathematical structure of the problem. 

To our knowledge, the only authors within the fi- 
nancial economics literature that have considered the 
possibility of no n-constant mortalit y rates in a life- 



cycle model are ICocco and Gomes (2009), although 



their Lee-Carter mortality model is not quite stochas- 
tic as in Milevskv and Promislowl (l200ll) . Dahl (2003). 



Cairns et al. or the various models described in 

the the book bv Pitacco et alj ( 2008 ). or the concerns 
expressed bv Norb erg (2010l). 

Moreover, a numbe r of very recent p aper s - for exam- 
ple iMenoncinI (l2008h . IStevens' (l2009h and iPostI (1201 oh 
- have examined the implications of (truly) stochastic 
mortality rates on the demand and pricing of certainly 
annuity products, but have not derived the impact of 
stochasticity on optimal consumption alone or exam- 
ined the impact of pure uncertainty in the mortality rate. 



Ano ther related paper is iBommier and Villeneuve 
(120 121) who examine the impact of relaxing the assump- 
tion of additively separable utility and what-they-call 
risk neutrality with respect to life duration. But, they 
also assume a deterministic force of moratlity in their 
formulation and examples. In that sense, our work is 
similar because we also relax the so-called risk neutral- 
ity and the intertemportal additivity. 

In sum, to our knowledge, none of the existing pa- 
pers within the LCM literature have assumed a stochas- 
tic force of mortality - which is the model of choice in 
the current actuarial and insurance literature - and then 
derived its impact on pure consumption behavior We 
believe this to be a foundational question, and in this 
paper our objective is straightforward, namely, to com- 
pare the impact of stochastic vs. deterministic mortality 
rates on the optimal consumption rate. 

1.1. A Proper Comparison 

Assume that two hypothetical retirees - i.e. con- 
sumers who are not expecting any future labour in- 
come - approach a financial economist for guidance on 
how they should spend their accumulated financial cap- 
ital over their remaining lifetime; a time horizon they 
both acknowledge is stochastic. Assume both retirees 
have time-separable and rational preferences and seek 
to maximize discounted utility of lifetime consump- 
tion under the same elasticity of intertemporal substi- 
tution (1/y), the same subjective discount rate (p) and 
the same initial financial capital constraint (Fq). They 
have no declared bequest motives and - for whatever 
reason - neither are willing (or able) to invest in any- 
thing other than a risk-free asset with instantaneous re- 
turn (r); which means they are not looking for guidance 
on asset allocation or annuities^] All they want is an op- 
timal consumption plan (c*(f); t >Q) guiding them from 
time zero (retirement) to the last possible time date of 
death (f < D). Most importantly, both retirees agree 
they share the same probability-of-survival curve de- 
noted by p{s). In other words they currently live in the 
same health state and the same effective biological age. 
For example, they both agree on a p(35) = 5% prob- 
ability that either of them survive for 35 years and a 
p{2Q) = 50% probability that either of them survive for 
2 years , etc. 



Yaaril (11964 Il965h showed exactly how to solve such 
a problem. He derived the Euler-Lagrange equation for 
the optimal trajectory of wealth and the related con- 
sumption function. 
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In Yaari's model both of the above-mentioned retirees 
would be told to follow identical consumption paths 
until their random date of death. In fact, they would 
both be guided to optimally consume c(f)* - F{t)/a{t), 
where a(t) is a function of time only and is related to an 
actuarial annuity factor We will explain this factor in 
more detail, later in the paper 

But here is the impetus for our comparison. Although 
both retirees appear to have the same longevity risk as- 
sessment and agree-on the survival probability curve 
p{s), they have differing views about the volatility of 
their health as proxied by a mortality rate volatility. In 
the language of current actuarial science, the first re- 
tiree (1) believes that his instantaneous force of mortal- 
ity (denoted by A^^^{t)) will grow at a deterministic rate 
until he eventually dies, while the second retiree (2) be- 
lieves that her force of mortality (denoted by A^^{t)) 
will grow at stochastic (but measurable) rate until a ran- 
dom date of death. As such, the remaining lifetime ran- 
dom variable for retiree 2 is doubly stochastic. While 
this distinction might sound farfetched and artificial, a 
growing number of researchers in the actuarial literature 
are moving to such models 0, rather than the simplis- 
tic mortality models traditionally used by economists. 
The actuaries' motivation in advocating for a stochastic 
force of mortality, is to generate more robust pricing and 
reserving for mortality-contingent claims. These studies 
have all argued that SfM models better reflect the uncer- 
tainty inherent in demographic projections vis a vis the 
inability of insurance companies to diversify mortality 
risk entirely. We ask: how do the recent actuarial mod- 
els impact the individual economics of the problem? 

When one thinks about it, real-life mortality rates 
are indeed stochastic, capturing (unexpected) improve- 
ments in medical treatment, or (unexpected) epidemics, 
or even (unexpected) changes to the health status of an 
individual. Rational consumers choosing to make sav- 
ing and consumption decisions using models based on 
deterministic mortality rates would likely agree to re- 
evaluate those decisions if their views about the values 
of those mortality rates change dramatically. Our thesis 
is that economic decision-making can only be improved 
if mortality models reflect the realistic evolution of mor- 
tality rates. 

We will carefully explain the mathematical distinc- 
tion between deterministic and stochastic forces of mor- 
tality (SfM) in Section |2] of this paper, but just to 
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make clear here, at time zero both our hypothetical re- 
tirees agree on the initial survival probability curve p(s). 
However, at any future time their perceived survival 
probability curves will deviate from each other depend- 
ing on the realization of the mortality rate between now 
and then. 

Motivated by such models of mortality, in this paper 
we derive the optimal consumption function for both re- 
tirees; one who believes in - and operates under a - 
stochastic mortality and one who does not. Stated differ- 
ently, we will solve the (consumption onlv) lYaari ( 1965 ) 
model where the optimal consumption plan is given as a 
function of wealth, time and the evolving mortality rate 
as a state variable. Indeed, with thousands of LCM pa- 
pers in the economic literature over the last 50 years, 
and the growing interest in stochastic mortality models 
in the actuarial community, we believe these results will 
be of interest to both communities of researchers. 

Recall that in the Yaari model conditioning on the 
mortality rate was redundant or unnecessary since its 
evolution over time was deterministic. All one needed 
was the value of wealth F{t) and time f. But, in a 
stochastic mortality model, the mortality rate itself be- 
comes a state variable. In this paper we show how 
the uncertainty of mortality interacts with longevity risk 
aversion (y) - which is the reciprocal of the intertempo- 
ral elasticity of substitution - to yield an optimal con- 
sumption plan. Mortality no longer functions as just a 
discount rate. 

To briefly preview our results, we describe the condi- 
tions under which retiree 1 (deterministic mortality) will 
start-off consuming more than retiree 2 (stochastic mor- 
tality), as well the conditions under which retiree 1 con- 
sumes less than retiree 2, and the (surprising) conditions 
under which they both consume exactly the same. We 
provide numerical examples under a variety of specific 
mortality models and examine the magnitude of this ef- 
fect. 

The remainder of this paper is organized as follows. 
In Section |2] we explain in more detail exactly how a 
stochastic model of mortality differs from the more tra- 
ditional (and widely used in economics) deterministic 
force of mortality. In Section [3] we take the oppo rtunity 
to review the (consumption only) lYaaril (119651) model 
and set our notation and benchmark for the stochastic 
model. In Section |4] we characterize the optimal con- 
sumption function in the stochastic mortality model un- 
der the most general assumptions, and prove a theorem 
regarding the relationship between consumption in the 
two models. In Section |5] we make some specific as- 
sumptions regarding the stochastic mortality rate and il- 
lustrate the magnitude of this effect, and Section |6]sum- 
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marizes our main results and concludes the paper The 
appendix contains mathematical details and algorithms 
that are not central to our main economic contributions. 

First, we explain exactly the difference between de- 
terministic and stochastic force of mortality. 



2. Understanding The Force of Mortality 

In most of the relevant papers in the LCM literature 
over the last 45 years the force of mortality from time 
zero to the last possible date of death is known with cer- 
tainty. Ergo, the conditional survival probabilities over 
the entire retirement horizon are known (in advance) at 
time zero. So, if a 65-year-old retiree is told (by his 
doctor) that he faces a 5% chance of surviving to age 
100 and a 37% chance of surviving to age 90, then 
by definition there is a 13.5% = (0.05/0.37) probabil- 
ity of surviving to age 100, if he is still alive at age 90. 
In other words, he makes consumption decisions today 
that trade-off utility in different states of nature, know- 
ing that if-and-when he reaches the age of 90, there will 
only be a 13.5% chance he will survive to age 100. In 
the language of actuarial science, the table of individual 
[q^+i; i - 0, . . . ,N] mortality rates is known in advance. 
This is the essence of a deterministic force of mortality 
and textbook life contingencies. If ^65 is the retiree's 
probability of dying between age 65 and 66, while ^66 
is the probability of the same retiree dying between age 
66 and 67, then the probability of surviving from age 65 
to age 67 is (1 - qeiXl - qee)- 

In stark contrast, under a stochastic force of mortal- 
ity the above multiplicative relationship breaks down. 
We do not know in advance how survival probabilities 
will evolve. While a 65 -year-old might currently face a 
5% estimated probability of surviving to age 100 and a 
37% chance of reaching age 90, there is absolutely no 
guarantee that the conditional survival probability from 
any future age, to age 100 (given the observed mortality 
rates), will satisfy the ratio. At time zero there is an ex- 
pectation of what the probability will be at age 90. But, 
the probability itself is random. This way of thinking - 
which might be new to economists - is the essence of a 
stochastic force of mortality and is the impetus for our 
paper. 

Here it is formally. Let A{t) denote the mortality rate 
of a cohort of a population, which may be stochastic or 
deterministic. Let f, = cr{A{q) | ^ < f) be the filtra- 
tion determined by A. Then individuals in the popula- 
tion have lifetimes of length ^ satisfying 



Assume further that A(t) is a Markov process, and define 
the survival function p(t, s. A) by 



pit, s,A)^E [e-f-^W''* I A(t) = a] . 



(2) 



This gives the conditional probability of surviving from 
time f to time s, given knowledge of the mortality rate 
at time f . Therefore 



> s \ ^ > t,r,) 



= £[e-r-^W''«ir,] = Mf,^,^(f))- (3) 

If f = then we write p(s, A) for p(0, s. A). 

Our basic problem in this paper will be to compare 
optimal consumption under two models that share a 
common initial value Aq of the mortality rate, as well 
as a common survival function p{t, Aq). Typically one 
will be deterministic and one stochastic. When we do 
actual computations, we will either choose a specific de- 
terministic model and calibrate a stochastic model to it, 
or conversely, we will choose a stochastic model and 
calibrate the deterministic model to it. Both possibil- 
ities are discussed below. It should be clear from the 
context which model we are discussing. But when it 
is necessary to make this distinction explicitly, we will 
write A°^{t) and A^^{t). 

2.1. Deterministic force of Mortality (DfM) 

Let Aq - A{<d) be the initial value of the mortality rate. 
In the deterministic case. 



p(f,^o) = e-/o■^('^)''«, 



(4) 



- A(q)dq 



(1) 



and we can recover A{t) as -pt(t, Ao)/p(t, Aq), where the 
f-subscript denotes the time derivative. In other words, 
if we start with a concrete stochastic model, and obtain 
the survival curve /?(?, /Iq) from it, the above formula 
determines the calibration of the deterministic force of 
mortality model. This approach is computationally sim- 
pler, but has the disadvantage that neither the stochastic 
nor deterministic model is in a simple form, familiar to 
and used by practitioners. In other words, a "simple" 
model for the stochastic force of mortality rates leads 
to a "complicated" model for the deterministic force of 
mortality, and vice versa. 

When doing actual calculations we will start by as- 
suming that A{t) follows a standard Gompertz model. 
The Gompertz model was introduced i n 1825, bu t 
more recently was popularized by Ca rrierel (Il994b . 
for example. Alternative models are presented in 
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Gutterman and VanderhoofI ( 19981) a nd others are dis- 
cussed as early as BrilUngen i 1961 ). In our case, we 



dA(t) - rjAit) dt 



(5) 



so A(t) = Aoe'''. The usual form for Gompertz is A(t) = 
l^-igix+i-m)/!}^ SO here we are using rj - Ijb and Aq - 
. This model is simple, and takes advantage 



-\g(x-m)lb 



of long experience calibrating the Gompertz model to 
real populations. 

Note that in the deterministic setting. 



p{t, s,A(t)) = g-f 

= e-/„'^(?)^?/e-/„'-i(?)rf9 = pCMo) 
piUAo)' 



(6) 



This will typically NOT be true in the stochastic setting. 
As long as we keep in mind that we are calibrating at 
time (i.e. to pit, A(,) only) that should not cause prob- 
lems. 

Table 1 displays a typical (loosely based on U.S. uni- 
sex annuitant mortality) deterministic mortality survival 
probability "matrix" of values together with the corre- 
sponding mortality rate at each age x, on the bottom 
row. Note that these numbers were generated using a 
(deterministic) Gompertz model in which m - 89.335 
and b - 9.5. Indeed, given the initial probability of sur- 
vival from age 65 to any age y > 65 (which is the first 
column in Table 1) one can solve for the conditional 
survival probability from age y to any age z > y, by di- 
viding the two probability values. This is the essence of 
equation Alas, when mortality rates are stochastic 
all numbers pit, s, /1(f)) beyond the first column in Table 
1, are unknown at time zero. 

2.2. Stochastic Force of Mortality ( SfM) 

There are many possible stochastic m odels to choose 
from. Starting from the models of Lee an d Carter 



1992h.lCaims et al.1 (l2006l) as well as lWills and Sheiiis 



20101) . actuaries have employed a variety of specifica- 



tions for the stochastic Ait), subsequently used to price 
mortality and longevity risk. In what follows in the nu- 
merical examples, we adopt a lognormal mortality rate, 
which is often called the Dothan model for interest rate s 
in the derivative pricing literature - see Dothanl (119781) . 



Although it might seem natural to have constant drift 
and diffusion coefficients, in order to calibrate to a given 
deterministic model, we allow a time-dependent growth 
coefficient. For most of the numerical examples pro- 
vided later-on we take: 



where Bit) is a Brownian motion. This is obviously the 
source of randomness in the stochastic force of mor- 
tality. There are many ways to select (or calibrate) 
a stochastic force of mortality to a particular survival 
curve. The details on how to actually compute this are 
provided in the second part of the appendix. 

With the probability background out of the way, we 
now review the (consumption only) |^ari||l965) model 
which is based on a deterministic force of mortality. 



3. Review of the lYaaril Il96^ Model 



The canonical lifecycle model (LCM) with a random 
date of death and assuming no bequest motive, can be 
written as follows: 



7 = max£'[ J" e '"uicit))l[,<^)dt^. 



(8) 



where ^ is the remaining lifetime satisfying Pr[^ > t\ — 
pit, A^)), defined above in Section|2l We fix a (determin- 
istic) last possible time D of death, so ^ < D. When 
the mortality rate is deterministic one can obviously as- 
sume independence between the optimal consumption 
c*(f) and the lifetime indicator variable l{t<Q, so that by 
Fubini's theorem we can re-write the value function as: 

J - max I e '"M(c(f))£'[l|,<^|]t/f 
Jo 

= max I e ^'ui 
" Jo 



'iicit))pit, Ao)dt. 



(9) 



From this perspective, there really is not any more ran- 
domness in the model. This is a problem within the 
calculus of variations subject to some constraints on the 
function c(f). In the end, the survival probability is ab- 
sorbed into the discount rate. 

Let r denote the risk free interest rate. To avoid 
the distractions of inflation models and assumptions, 
throughout this paper we assume that r is expressed in 
real (after-inflation) terms and therefore consumption 
c(f) is expressed in real terms as well. The wealth (bud- 
get) constraint can then be written as: 



Frit) = rFit) + no- c(f), 



(10) 



dAit) - fiit)Ait) dt + o-Ait) dBit) 



(7) 



with boundary conditions F(0) = W > and FiD) - 0. 
We are using the subscript F, to denote a first derivative 
w.rt time, and if needed F„ for the second derivative. 
The parameter denotes a constant income rate which 
we include in this section for comparison with Yaari's 
model, but which in subsequent sections will be taken 
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Table 1: Conditional Survival Probability: Deterministic Mortality 




x = 65 


X = 70 


X = 75 


X = 80 


X = 85 


X = 90 


X = 95 


X = 100 


To Age 65 


1.000 
















To Age 70 


0.9479 


1.000 














To Age 75 


0.8659 


0.9135 


1.000 












To Age 80 


0.7429 


0.7837 


0.8580 


1.000 










To Age 85 


0.5733 


0.6047 


0.6620 


0.7716 


1.000 








To Age 90 


0.3696 


0.3899 


0.4268 


0.4975 


0.6447 


1.000 






To Age 95 


0.1758 


0.1855 


0.2031 


0.2367 


0.3067 


0.4757 


1.000 




To Age 100 


0.0500 


0.0527 


0.0577 


0.0673 


0.0872 


0.1353 


0.2844 


1.000 


A(x) 


0.0081 


0.0137 


0.0232 


0.0394 


0.0667 


0.1129 


0.1911 


0.3234 



to equal zero; c{t) is the consumption rate and the con- 
trol variable in our problem. In a follow-up paper we 
hope to examine the impact of additional factors, such 
as different interest rates for borrowing versus lending, 
of the availability of actuarial notes (i.e. the case when 
the interest rate is r + A(t)). 

In this paper we operate under a constant relative risk 
aversion (CRRA) formulation for the utility function. In 
principle this should mean using u(c), where: 



u(c) - 



-I 



1 



■7 



(11) 



for 7 > and y + \, with the understanding that when 
y = 1 we define m(c) = Inc. This family of utilities 
varies continuously with y. The marginal utility of con- 
sumption is the derivative of utility with respect to c, 
which is simply 



> 0. 



(12) 



Of course, it makes no difference to our optimization 
problem (and the optimal control) if we shift u by an 
arbitrary additive constant. So to make scaling relation- 
ships easier to express, actual calculations will be car- 
ried out using the equivalent utilities 



m(c) - 



1-r 



(13) 



for y > and j + \. When y = 1 we take m(c) - u(c) - 
Inc. 

As a consequence of the Euler-Lagrange Theorem, 
the optimal financial capital trajectory F{t) must satisfy 
the following hnear second-order non-homogenous dif- 



ferential equation over the values for which F{t) + 0. 

f „(„ - (l;£^ . ,) f ,,„ + .(l^i^) f (,) 

IZEziW),.. ,,4) 

When the pension income rate ;7ro = the differen- 
tial equ a tion c ollapses to the homogenous case. See 
Chian^l (1992), for examplem for an exposition of 



Euler-Lagrange equations in economics. 

3.1. Explicit Solution: Gompertz Mortality 

When the (deterministic) mortality rate function 
obeys the (pure) Gompertz law of mortality 

1 I X + t — m\ 
'^(0 = ^exp(— (15) 

the survival probabihty can be expressed as 

p{t,XQ) = e-/o*^)«'? = gWod-^'"")^ (16) 

Here x denotes the age at time 0, m is called the modal 
value and h is the dispersion coefficient for the Gom- 
pertz model. To simplify notation let 



k(t) 



r-p-A(t) 
7 



(17) 



and recall from the budget constraint that: 

c(f) = rFit) - F,{t) + JTo, 
c,(f) = rF,(t) - F„{t). 

Equation (fT4l i can be rearranged as 

F„{t) - rF,(t) + k(t)(rF(t) - F,(t)) = -^(0^0, (20) 



(18) 
(19) 
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which then leads to 

kityit) - c;{t) = 

The solution to this basic equation is 



(21) 



= c*(0)eT^'e"r /o'-'W''-' = c{0)e^'p{t, Aq)^''', (22) 

where c*(0) is the optimal initial consumption rate, to 
be determined, which is the one free constant resulting 
from equation (1211 1. Note that when the interest rate r 
is equal to the subjective discount rate p, and y - I 
(i.e. log utility), the optimal consumption rate at any 
age X + f is the probability of survival to that age times 
the initial consumption c*(0). However, when y > 1, 
which implies higher levels of risk aversion, the opti- 
mal consumption rate will decline at a slower rate as the 
retiree ages. Longevity risk aversion induces people to 
behave as if they were going to live longer than deter- 
mined by the actuarial mortality rates. We will explore 
the impact of y on the optimal consumption path in a 
stochastic force of mortality model, later in Section 2] 
which is why it's important to focus on this here. 



Mathematically one can see that (p(t,Ao)) 



I/(y+£) 



IS 



greater than {p(t, /Lo))'^'' for any e > since p{t, Aq) < 1 
for all t. Finally, note that in the Gompertz mortality 
model evaluating (p(t, /lo))'^'' for a given (x, m, b) triplet 
is equivalent to evaluating pit, Aq) under the same x, b 
values, but assuming that m* - m + blny. This then 
implies that one can tilt/define a new deterministic mor- 
tality rate Ao - yAo and derive the optimal consumption 
as if the individual was risk neutral. This will be used 
later in the explicit expression for F(t) and c*{t) 

Moving on to a solution for F(t), we now substi- 
tute the optimal consumption solution ( l22l i into equa- 
tion (fTST l to arrive at yet another first-order ODE, but 
this time for F{ty. 

F,(t) - rF(t) -710 + c*(0)e'T'p(t, Ao)"^ = 0. (23) 

Writing down the canonical solution to this equation 
leads to: 



^ e-T'p{s, Ao)^'^e-'''ds + F{0) | , (24) 



Fit) = 

-c*(0) 



where Fq denotes the free initial condition from the 
ODE for Fit) in equation ( l23l l. Recall that we still 
haven't specified c*(0), the initial consumption). We 



will do so (eventually) by using the terminal condition 
FiD) = 0. 

To represent the wealth trajectory explicitly define the 
following (new) Gompertz Present Value (GPV) func- 
tion 



alir,m,b)^ pis, Ao)e-''ds ^ K^'^^^'^^'^'d. 
Jo Jo 



= r e *V*' r'ds (25) 
Jo 

bVi-rb, exp{^}) - bTi-rb, exp{^!^^}) 
exp{(m - x)r - exp{^)} 

The function a].ir, m, b) - ait) is the age-x cost of a 
life-contingent annuity that pays $1 per year continu- 
ously provided the annuitant is still alive, but only until 
time t - T, which corresponds to age x + T. If the in- 
dividual survives beyond age (x + T) the payout stops. 
Naturally, when T - co the expression collapses to a 
conventional single premium income annuity (SPIA). 

Note that r(A, B) is the incomplete Gamma function. 
In other words, equation (l25t is analytic and in closed- 
form. 

The reason for introducing the GPV is that combining 
equation (l24l l with equation (|25] | leads to the (very tame 
looking) expression 



Fit) = (piO) + ^ j e" - a'jr - k. 



m*,fo)c*(0)e''' - — , 
(26) 

where recall that m* - m+blny. Then, using the bound- 
ary condition Fr - 0, where r is the wealth depletion 
time, we obtain an explicit expression for the initial con- 
sumption 



c'iO) 



iFii)) + nolr)e'' -mlr 
a\ir - k, m*, b)e'"^ 



ill) 



3.2. Consumption Under DfM: Numerical Examples 

In our numerical examples we assume an 86.6% 
probability that a 65-year-old will survive to the age of 
75, a 57.3% probability of reaching 85, a 36.9% prob- 
ability of reaching 90, a 17.6% probability of reaching 
age 95 and a 5% probability of reaching 100. These 
are the values generated by the Gompertz law with 
m - 89.335 and b - 9.5. To complete the param- 
eter specifications required for our model, we assume 
the subjective discount rate (p) is equal to the risk-free 
rate r - 2.5% Within the context of a Ufecycle model, 
this implies that the optimal consumption rates would 
be constant over time in the absence of longevity and 
mortality uncertainty. 
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We are now ready for some results. Assume a 65- 
year-old with a (standardized) $100 nest egg. Initially 
we allow for no pension annuity income (tto = 0) and 
therefore all consumption must be sourced to the invest- 
ment portfolio which is earning a deterministic interest 
rate r = 2.5%. The financial capital F(t) must be de- 
pleted at the very end of the lifecycle, which is time 
D = (120 - 65) = 55 and there are no bequest motives. 
So, according to equation dZTl l. the optimal consumption 
rate at retirement age 65 is $4,605 when the risk aver- 
sion parameter is y = 4 and the optimal consumption 
rate is (higher) $4,121 when the risk aversion parameter 
is set to (higher) y = 8. 

As the retiree ages (f > 0) he/she rationally consumes 
less each year - in proportion to the survival probability 
adjusted for y. For example, in our baseline y = 4 level 
of risk aversion, the optimal consumption rate drops 
from $4,605 at age 65, to $4,544 at age 70 (which is 
t - 5), then $4,442 at age 75 (which is f = 10), then 
$3,591 at age 90 (which is t = 25) and $2,177 at age 
100 (which is f = 35), assuming the retiree is still alive. 
A lower real interest rate (r) leads to a reduced optimal 
consumption/spending rate. All of this can be sourced 
to equation (l22l l. 

Thus, one of the important insights is that a fully ra- 
tional consumer will actually spend less as they progress 
through retirement. The optimizer spends more at ear- 
lier ages and reduces spending with age, even if his/her 
subjective discount rate (SDR) is equal to (or less than) 
the real interest rate in the economy. 

Intuitively the individual deals with longevity risk by 
planning to reduce consumption - if that risk material- 
izes - in proportion to th e surv i val pr obability, linked 
to their risk aversion. The lYaaril(ll965l) model provides 
a rigorous foundation to the statement by iFisher (Il930h 
in his book Theory of Interest (page 85); ". . . The short- 
ness of life thus tends powerfully to increase the degree 
of impatience or rate of time preference beyond what 
it would otherwise be. . . " and (page 90) "Everyone at 
some time in his life doubtless changes his degree of im- 
patience for income. . . When he gets a little older,. . . he 
expects to die and he thinks: instead of piling up for the 
remote future, why shouldn't I enjoy myself during the 
few years that remain ? '0 



'For additional (case specific) examples of the lYaaril i 19651) model 
in action during the non-labour income retireme nt phase, we refer the 
interested reader to iMilevskv and Huand | |2010|) or a recent paper by 
iLachancd l20ll) . 



3.3. Time-zero Consumption Ratio — Initial Withdrawal 
Rate 

Finally, in the very specific case when ttq = (which 
implies that the wealth depletion time is t - D) and 
the subjective discount rate p - r, the retiree must rely 
exhaustively on his/her initial wealth Fq. We get 



1 



c*(0) ^ 

F(Q) aD(r-AQ/y,m*,b) 



(28) 



We now have all the ingredients to compare with a 
stochastic model. This ratio is often called the Initial 
Withdrawal Rate (IWR) amongst financial practitioners 
and in the retirement spending literature. 

4. Optimal Consumption: General Results 

In this section we obtain the most general optimal 
consumption strategy for a retiree maximizing expected 
discounted utility of consumption with uncertai n life- 
time, w hich will include the (consumption only) lYaari 
(Il965h model as a special case. Since our main focus 
now is on the mortality model, at this stage we make 
the additional assumption p - r, that is, that the subjec- 
tive discount rate equals the interest rate in the economy. 
Also, in contrast to the discussion in the previous sec- 
tion, we assume no exogenous pension income, so that 
ttq - 0, which then precludes any borrowing. We now 
assume a fixed terminal horizon T, which denotes the 
last possible date of death. The mathematical formula- 
tion is to find 



J = max E 

c(s) adapted 



g- ];('■+■*(?)) ''9„(c(s))c/s 

Jo 



A{Q) = A, F(Q) = F 



(29) 



WTiereas in Section |3] of this paper we used calculus of 
variation techniques to derive the optimal trajectory of 
wealth and the consumption function, given the inclu- 
sion of mortality as a state variable we must resort to 
dynamic programing techniques to obtain the optimaUty 
conditions. Regardless of the different techniques, we 
will show h ow the optim al consumption function col- 
lapses to the lYaari d 19651) model when the volatility of 
mortality is zero. 
Define: 



J(t, A, F) - max E 

c{s'} adapted 



XT 
g-J^ '''?M(c(i))t/i 



A(t) = A, Fit) = F 



(30) 
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As in the deterministic mortality model, the wealth pro- 
cess (which we shall soon see is stochastic) satisfies 
dF{t) = {rF{t) - c(t))dt. Assume that there is an op- 
timal control. Then for that control. 



IX' 



+ f e-.r(^+-^(''»'''?M(c(i))c/s (31) 
Jo 



is a martingale. This will likewise give a super- 
martingale under a general choice of c. Applying 
Ito's lemma, we obtain the following Hamilton- Jacobi- 
Bellman (HJB) equation: 



sup {u(c) - cJf] + J, -(r + A) J 



2)2 



+ rFJf + iJ(t)AJ,i + 



-Jaa^O. (32) 



If there is any possibility of confusion, we will denote 
this value function y^™(r, A, F). 

For deterministic mortality, HJB can be obtained by 
sendi ng cr — > w ith //(?) - rj, which was equal to l/b 
in the lYaaril(ll965i) model derived in Section |3] as 



sup{u(c) - cJF} + Jt-(r + A)J + rFJf + T]AJA =0. (33) 

c 

Moving on to the optimal consumption plan, we solve 
the HJB equation under CRRA utility as follows: let 



u(c) 



1-r 



J 



a(t. A), 



(34) 



where the second expression results from the scaling 
which follows from the first, and apply the first order 

„ ' . We obtain c* - Fa r and get the 



condition c* 



following equation for a{t. A): 



1 . (T^A^ 

— (ry + A)a +ya y + ^{t)AaA H — a,u - (35) 



with boundary condition a{T, A) - 0. 

We now solve the PDE for a{t. A), which we re-write 
as: 

li, + l-L+-)l5+tx(t)AliA + ^^o-''A^I3l + ]-(T^A^Pu = 0. 
^ y' Ip I 

(36) 

for ;0 = /3{t,A) - a(t,Ay^'>' . The boundary conditions 
are (^(T, A) = 0, /3i(t, oo) = and at /I = we solve 
P, + I - rfi - Q. Note that the optimal consumption rate 
is c = Fip, using shorthand notation. On to the main 
theorem. 



4.1. Stochastic Force of Mortality: Main Theorem 

Denote by c^^{t. A, F) the optimal consumption at 
time f, given A(t) - A and Fit) - F, under a stochastic 
force of mortality (SfM) model. Denote by c°^(t, F) 
the optimal consumption at time t, when F{t) = F, un- 
der a deterministic force of mortality (DfM) model. 

Theorem 1. Assume that the survival functions for the 
two models agree: p^-^(t, Aq) = p'^^it, Aq) for every 
f > 0, and that utility is CRRA{y) 

(a) y > 1 ^ c^fl^iO, Aq, F) > c^^(0, F); 

(b) y = 1 ^ c^^^(0,Ao,F) = c^^iO, Fl- 
ic) < r < 1 ^ c^^iO, Ao, F) < c^-^(0, F). 

Proof. To see this, we change point of view, and work 
exclusively with the stochastic model. So we drop 
the SfM superscript, and write p - p^^, J - 7^™, 
c* - c^'^, A - /i^''^, etc. Within that model, we pose 
two different optimization problems, depending on the 
level of information available about Ait). The value 
function Jit, A, F) solves the problem given before in 
( |30] |, where c(f) can be any suitable process adapted to 
But we define a new value function 7'(f, F) in which 
we impose an additional constraint on c(0, namely that 
it be deterministic. More precisely. 



7(0, Ao, Fq) - max E 

c{s) adapted 



and 



nT 

e-.fo('-^'*(i^1diuicis))ds 

Jo 



(37) 



j\0,Fo) (38) 

r-T 

Jo 

e^'^^pis, Ao)uicis))ds. 



max E 

cis) deterministic 



max 

c{s) detenninistic 



Jo 



We let c* denote the optimal control for J, and c' denote 
the optimal control for 7 ' . 

Since every deterministic control c(f) is also adapted, 
we have the basic relationship 



JiO,Ao,Fo)>riO,Fo). 



(39) 



On the other hand, the above expression is exactly what 
the old deterministic model would have given. That is. 



7'(0,Fo) = 7°^(0,Fo) 



(40) 



„DfM 



and c' - c 

Due to scaling, Jit,A,F) = ait, A)F^-^ /il - y) and 
c*(f. A, F) - ait, A)^^^'>'F for some function a > 0. Like- 



wise J^^it,F) = aiit)F^-^/il -y) and c' 



-i/r 



F for 
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someoi > 0. Ify > 1 then 1-y < 0, sofl(0,/l()) < fli(O), 
so c* > c' at f = 0. This shows (a). Likewise if 
< 7 < 1 then a(0,Ao) > fli(O), so c* < at f = 0. 
This shows (c). 

Recall that when -y - 1, we have u(c) - Inc. Earlier, 
when 7 9^ 1, we had u(c) - c'"''/(l - j) and could 
make use of a scaling relation. In other words, if c is 
optimal for F, then kc is optimal for kF, and that leads 
to the expression J(t, A, kF) = k^^~^J{t, A, F). Or in other 
words, 

J{t,A,F)^F^-^J(t,A,\). (41) 

With logarithmic utility, the corresponding expression 
is that 

J{t, A, kF) = J(t, A, F) + (In k) f e-'^'-'^pit, s. A) ds. Or 
in other words, 

J{t,A,F)^ J(t,A,l) + (In F) J e-'^'-''> p(t,s,A)ds. 

(42) 



Likewise, 



7°^(f,F) 



7°™(f, 1) + (1 



pit,Ao) 



ds. 
(43) 



The first order conditions in the optimization problem 
then imply that 



c* = F/ 



I- 



-r(s-t) 



p{t, s. A) ds. 



c^^^FI f e-'-^-^-'^P^ds. 
J, pit,Ao) 

These agree when we send f — > 0, showing (b). 



(44) 
□ 



The theorem certainly proves that y = 1 is a point of 
indifference. The invariance of mortality volatility when 
utility is logarithmic is reminiscent of similar results in 
consumption theory where income negates substitution 
effects. More on this later 

Note that we only use j\0, F) above, not y'(f, F). If 
we had, we would have had to be careful. The correct 
definition is that 



J'(t,F)^J°^(t,F) 



= max r e-'-'-'^-'^E^^u(c(s))ds (45) 
c(0 Jt p(t,Ao) 



rather than 



maxii 



J g-f ('•+'!(?» rf?„(c(5))t/i 

= max e-''''-'^E[p(t, s, A(t))]u(c(s))ds. (46) 



These quantities have connections to annuities, as sug- 
gested by the fact that the optimal consumption rates 
given above are, as a fraction of wealth, inverse an- 
nuity prices. In particular, e ''^p(s,Ao)ds is the 
(actuarial) price of a deferred annuity, purchased at 
time with payments starting at time t. While 
e "£■[/?(?, s, A(t))] ds is a forward annuity price. 
That is, if at time an insurance company guarantees 
(a retiree) the right to buy an annuity at time f at a price 
determined at time 0, then this is that price (computed 
actuarially, i.e. by discounting mean cash flows). 

4.2. Intuition and Relation to Known Results 

How should one interpret our result? It is tempting 
to view stochastic mortality as simply "more risky" that 
deterministic mortality, but that is not in fact the rea- 
son consumption shifts. The true explanation for our 
result is that the comparison can be reinterpreted equiv- 
alently as one between two different control problems, 
both within the context of the stochastic hazard rate 
model. Namely, a control problem in which the haz- 
ard rate is observed, so one can react to changes, versus 
one in which the hazard rate is not observed, so the con- 
trol must be determined in advance. The utility in the 
deterministic model is the same as the utility for the sec- 
ond control problem (and indeed, this is the basis of our 
proof). So the mere presence of stochastic hazard rates 
will not cause a change in consumption; what shifts con- 
sumption is the ability to react to those changes. 

There are two possible reactions to that ability to ad- 
just consumption. One is to shift consumption into the 
future, taking advantage of the ability to adjust con- 
sumption upwards later, if the hazard rate should climb 
more than expected. The other reaction is to opt to con- 
sume more now, in the knowledge that one can cut back 
later if it seems likely that one will live longer than ex- 
pected. Our message is that either reaction can be ra- 
tional, and that which one is adopted depends on the 
person's risk aversion, with the switch occurring at the 
point of logarithmic utility. The choice is between act- 
ing more conservatively in view of the possibility one 
might live longer, versus acting more aggressively in 
view of one's ability to react to changes in the hazard 
rate. 

There are other results in the literature where loga- 
rithmic utility is a qualitative point of indifference in 
behavior. An example - in a completely different con- 
text - is the classical result on the equilibrium pricing of 
assets derived by .Lucas. (,1978,F i In a Lucas-type model 



We thank Thomas Davidoff for pointing out this analogy 
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- under logarithmic utility preferences - the equilibrium 
price of trees (or any other income producing asset for 
that matter) does not depend on the projected level of 
fruit output from those trees. The economic reason for 
this is that there are two effects on the current equilib- 
rium price, of an increase in the expected future amount 
of fruit from trees. The first is the fact that at any given 
marginal utility of consumption of the fruit, the higher 
expected fruit production increases the attractiveness of 
owning trees today, which raises the current price of a 
tree. But, at the same time, the increased expected fruit 
output in future periods means higher consumption and 
lower marginal utility of consumption in that future pe- 
riod. This effect tends to reduce the attractiveness of 
owning trees today - the tree is going to pay off more in 
a time when marginal utility is expected to be low - and 
thus lowers the current value (and hence price) of a tree. 
These two forces are the manifestation of the (pure) in- 
come effect and substitution effect from the theory of 
consumer choice, and their net result - i.e which actual 
dominates - depends on the shape (and curvature) of the 
utility function. 

In the case of logarithmic utility, income and substi- 
tution effects are of the same size and opposite sign so 
the two forces exactly offset each other, leaving the cur- 
rent price of a tree unchanged in the face of a rise in 
expected f uture fruit ou tput. This is (one of) the re- 
sults from iLucasI (119781) . Although it does not appear 
the same powers are at force in our stochastic mortality 
model, this does illustrate that there are a number of set- 
tings in which one finds that logarithmic utility (y = 1) 
is the point of indifference between two opposing con- 
sumption effects. 

5. Optimal Consumption: Numerical Examples 

We started with a particular survival probability at 
time zero, namely the Gompertz mortality curve with 
parameters m - 89.335 and b = 9.5. The age x = 65 
survival probabilities to any age y > x are given in Ta- 
ble 1 . Both hypothetical retirees agree on these num- 
bers, which means that their initial mortality rate is 
^0 = (1/9.5) exp{(65 - 89.335)/9.5) = 0.008125. 

Over time retiree 1 believes his mortality rate will 
grow at a rate rj - (1/9.5) - 0.105 263 16 per year, 
while retiree 2 believes it will evolve stochastically with 
a time-dependent growth rate of fi(t) and a volatility 



cr. The actual curve yu(f) depends on the selected pa- 
rameter for volatility, since fi{t) is constrained to match 
p{0, Aq). The actual process for extracting //(f) for any 
given value of cr is rather complicated (although it is not 
central to our analysis) and is placed in the appendix of 
this paper. With these numbers in hand - and specifi- 
cally the function //(f) for the drift of the mortality rate - 
we can proceed to solve the PDEs given in equation ( l35T l 
and (|36] |. which then lead to the desired optimal con- 
sumption function and the initial portfolio withdrawal 
rate at age 65. 

Table 2 provides a variety of numerical examples 
across different values of (mortality volatility) cr and 
(risk aversion) y, once again assuming that the retirees 
are both at age x = 65 with observable mortality rate 
Ao - 0.0081. As we proved in Section[3] and discussed 
above, the consumption rate is the same across all levels 
of mortality volatility when y - 1 . It increases relative 
to DfM when y > 1 and decreases relative to DfM when 
y < I. Notice the impact of stochastic mortality on op- 
timal withdrawal rates is reduced as the value of risk 
aversion increases. Notice how at a coefficient of rela- 
tive risk aversion y = 10, the portfolio withdrawal rates 
are approximately 4.6% at all listed volatility levels. 

Note that the cr values provided are rather ad hoc 
and have not been estimated from any particular de- 
mographic dataset. We refer th e interested r eader to 
recent actuarial papers - such as iBauer et al.l (12008 1 - 



to the authors, and we refer the interested reader to the lec- 
ture notes by Christopher Caroll, available at the following web- 
site: econ.jhu.edu/people/ccarroll/public/lecturenotes. 
for this particular interpretation (and derivation) of the Lucas model. 



for an empirical discussion around the estimation of the 
volatility of mortality. Our objective here is to explore 
whether or not mortality volatility has a (noticeable) im- 
pact on rational behavior as opposed to on insurance 
pricing. 

To sum up, when the coefficient of CRRA (denoted 
by y) is equal to one, and the retiree has logarithmic 
utility preferences, the optimal consumption rate at time 
zero is identical in both models. In other words, a retiree 
who cannot adjust their consumption plan as mortality 
rates evolve starts-off with the exact same consump- 
tion rate as the (more knowledgeable) consumer who 
can adapt to changes in mortality rates and health sta- 
tus. Although the path of their respective consumption 
will diverge over time - depending on the evolution of 
mortality rates - initially they are the same. In contrast, 
when the coefficient of CRRA is greater than one and 
the retiree is more risk-averse compared to a logarithmic 
utility maximizer, the initial consumption rate is higher 
in the stochastic model vs. the deterministic model. In 
other words, as one might expect the ability to adapt 
to changes in health status and new information about 
mortality rates allows the retiree to be more generous at 
time zero. Finally, when the coefficient of the CRRA is 
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Table 2: Optimal Retirement Portfolio Withdrawal Rates c*(0)/F() 


Mortality Volatility 


r = 0.5 


r = 1.0 


r = 1.5 


r = 3 


r = 5 


r = 10 


cr = 


7.59% 


6.12% 


5.58% 


5.02% 


4.78% 


4.61% 


a- = 15% 


7.52% 


6.12% 


5.60% 


5.04% 


4.80% 


4.62% 


a- = 25% 


7.44% 


6.12% 


5.62% 


5.06% 


4.82% 


4.63% 


Notes: Retirement age 65, interest rate r - 2%, mortality Aq - 0.0081 



between one and zero, the result is reversed. The canon- 
ical retiree in a stochastic mortality model will consume 
less compared to their neighbor who is operating under 
deterministic mortality assumptions. 

Not withstanding the above results, the absolute con- 
sumption rate at time zero is uniformly higher the lower 
the coefficient of relative risk aversion. This is a mani- 
festation of longevity risk aversion. The retiree is con- 
cerned about living a long time, and therefore consumes 
less today to protect themselves and self-insure con- 
sumption in old age. 

6. Discussion and Conclusion 

In this article we extended the lifecycle model (LCM) 
of consumption over a random-length lifecycle, to a 
model in which individuals can adapt behavior to new 
information about mortality rates. The lifecycle model 
of saving and consumption continues to be very popu- 
lar as a foundation model for decison-making amongst 
financial advisors, as recen tly described in the mono- 
graph by'Bodie et al.'('2008'). 

[Yaari ( 1964. 1965) was the first to include lifetime 
uncertainty in a Ramsey-Modigliani lifecycle model 
and amongst other results, he provided a rigorous foun- 
dation for Irving Fisher's claim that lifetime uncertainty 
increases consumption impatience and is akin to higher 
subjective discount rates. When the mortality rate itself 
is stochastic, this analogy is no longer meaningful and 
- to our knowledge - the pure lifecycle model has not 
been extended into the realm of 21st century models of 
mortality and longevity risk. 

We buUt this extension by assuming that (i) the in- 
stantaneous force of mortality is stochastic and obeys 
a diffusion process as opposed to being deterministic, 
and (ii) that a utility-maximizing consumer can adapt 
their consumption strategy to new information about 
their mortality rate (a.k.a. current health status) as it 
becomes available. Our diffusion model for the stochas- 
tic force of mortality was quite general, but inspired by 
(a.k.a. borrowed from) the recent literature in actuar- 
ial science. We focused our modeling attention on the 



retirement income stage of the LCM where health con- 
siderations are likely to be more prevalent and to avoid 
compHcations induced by wages, labor and human cap- 
ital consideration. 

In the first part of this paper we re-derived the opti- 
mal consumption function under a deterministic force 
of mortality (DfM) using techniques from the calculus 
of variations. We provided a closed-form expression for 
the entire consumption rate function under a Gompertz 
mortality assumption. With those benchmark results in 
place, we derived the optimal consumption strategy un- 
der a stochastic force of mortality (SfM), by express- 
ing and solving the relevant Hamilton-Jacobi-Bellman 
(HJB) equation. In addition to the time variable, two 
state variables in the resulting PDE are current wealth 
and the current mortality rate. 

Retirees with (i) no bequest motives, (ii) constant rel- 
ative risk aversion (CRRA) preferences, and (iii) subjec- 
tive discount rates equal to the interest rate are expected 
to consume less as they age since they prefer to allocate 
consumption into states of nature where they are most 
likely to be alive. This is the conventional diminish- 
ing marginal utility argument. In our model, a positive 
shock to the mortality rate in the form of pleasant health 
news (perhaps a cure for cancer) will reduce consump- 
tion instantaneously and further than expected at time 
zero. A negative shock to the mortality rate (for exam- 
ple, being diagnosed with terminal cancer) will increase 
consumption beyond what was expected. 

Moving forward, a natural extension would be to ex- 
plore the impact of stochastic investment returns as well 
as mortality rates an d include a str ategic asset alloca- 
tion dimension, a la iMertonI (Il97lb . Another item on 
our research agenda is to explore the optimal allocation 
to health and mortality-contingent claims in a stochastic 
mortality mo del. Recall that one of the noted results of 
Yaaril(ll965h is that lifecycle consumers with no bequest 



motives should hold all of their wealth in actuarial notes. 
However, in the presence of a stochastic mortality, it is 
no longer clear how an insurance company would price 
pension annuities, given the systematic risk involved. In 
such a model, a retiree would have to choose between 



12 



investing wealth in a tontine pool, with corresponding 

stochastic returns or purchasing a pension annuity with 
a deterministic consumption flow, but possibly paying 
a risk-premium for the privilege. We conjecture that 
in a stochastic mortality framework, the optimal prod- 
uct allocation is a mixture of participating tontines and 
guaranteed annuities. 
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Appendix A. Matching Time-Zero Survival Curves 

The caUbration of our economic model leads to an 
interesting by-product problem in actuarial finance. In 
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particular, in order to construct a stochastic force of 
mortality that matches or fits a pre-determined Gom- 
pertz survival curve - the most popular and frequently 
used analytic law in this literature - one requires a log- 
normal diffusion process in which the drift itself grows 
even faster than exponentially over time. In this ap- 
pendix we explain the mechanics of the procedure. 

Given a deterministic model (Gompertz in our nu- 
merical examples), we compute the time-zero survival 
function p{t, io). We match this using a stochastic 
model, by a suitable choice of parameters. This means 
that at time the two models deliver identical survival 
probabilities. Recall that at times other than t = the 
comparison will no longer be meaningful, even con- 
trolling for the current observed mortality rate, because 
the mismatch between conditional survival probabilities 
means that the two models give different views of life- 
times going forward. 

Let A{t) = e"/o''<«^''? and define a pseudo-density 
q{t, A) by the formula 

(p(A)q(t,A)dA. (A.l) 



Then p{t, Ao) = q(t. A) dA. By Ito's lemma, 

^W(0)A(/) = 0(^) + ^ K{s)\^{s)A{s)(l)'{A{s)) 



+ —A{,sY^"iA{,s))-A{s),t>im) 



ds 



f 

Jo 



A{sMs)A{s)(p'{Ais))dB{s). (A.2) 



Take expectations and differentiate with respect to t. We 
get 



r cf>(A)q,(t,A)dA= f \tim<P'(A) 
Jo Jo 



+ —A'<l>"iA)-A<f>iA) 



qit,A)dA (A.3) 



with initial condition q(0, ■) = 6^. Using integration by 
parts (for vanishing fast at and oo), we have 



qt(.t,A) 



+ Y-g^[^^q(f'^-Mt,^)- (A.4) 



So if //(f) is known for < ? < ?i, then all expectations 
^ q(ti, A)(f)(A) dA can be found by solving the forward 
equation for q and then integrating against (f>. 



Let = AqdA and A(2) = £° A^qdA be the first 
two moments of q(t, A). Note that the zeroth moment is 
the survival probability, so we can integrate (by parts) 
the forward PDE for q and the product of A and the for- 
ward PDE and obtain the following relationships 



dp 

^^(2) = Kt)A^i) - 



dA: 



(1) 



dt 



Combined the two expressions, we have 



ix(i) = 



A(V) 



(A.5) 



(A.6) 



Replacing /i(f) in the forward PDE for q and obtain an 
integro-differential equation 



qt{t,A) = - 



A(i) dA 
a q2 



[Mit,A)] 



+ YdA^ [-I'^C?. ^)\ - Mit, (A.7) 



or 



qtit,A) 



-'^ + j^'°A^q(t,^d_ 
dA 



dp 

dt 



[Aq(t,A)] 



+ y ^ [A^qit, A)] - Aq{t, A), (A.8) 

which we can solve numerically with the initial condi- 
tion ^(0, A) = 6(A - Aq). 

We solve the integro-differential equation for q nu- 
merically first, obtain the value of ^{t). We then solve 
the HJB equation for optimal consumption as before, 
with the constant ju now replaced by the function fi{t) at 
A = 0. 

Finally, we should record a couple of remarks about 
the form of fi{t). First of all. 



m = v. 



(A.9) 



To see this, observe that pt{Q, /to) = -E^AqKq] = -Aq. 
So Al = E[Alhl\ = p„(Q,Ai{) + M,p(Q). But p„(0,.lo) 
can be computed exphcitly, since it is Gompertz, to give 
/Iq - Aqti. This impUes that //(O) = r]. 

Second, note that //(f) should be increasing in cr. The 
mean £[e-/o' does not change with cr, so by con- 
vexity of the exponential, the median of this quantity 
must decrease as we increase the variance. In other 
words, n{t) must rise. Put another way, this expecta- 
tion is driven by the possibility of relatively larger val- 
ues of the exponent, ie of abnormally low values of A. 
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As a rises, the impact of longevity risk gets more pro- 
nounced, and to compensate for that the growth rate 
must also rise. 
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